TECHNION-PHYS-94-16 

E^p-ph/ 941 1422 



Soft gluon suppression of l/N^ contributions in 
^ ; color suppressed heavy meson decays 

>: 

yyr '. Igor Halperin 

Q\ I Technion - Israel Institute of Technology 

^ ■ Department of Physics 

^ '. Haifa, 32000, Israel 

K^ . e-mail address: hieorOtechunix. technion. ac.il 

(N. 

^' 

^ : Abstract: 

"^ ! We discuss the non-factorizable terms in color suppressed (Class II) decays. Our 

CL|| emphasis is on the non-perturbative soft gluon exchange mechanism, which has 

p • been previously found to be responsible for the rule of dicarding 1/Nc in the Class 

^ ■ I decays. The non-factorizable contribution to the decay B^ -^ D^tc^ at the tree 

level is estimated within the light cone QCD sum rule method which combines the 

. ^ ! technique of the QCD sum rules with the description of the pion in terms of the set 

rS [ of wave functions of increasing twist. We find that the same soft gluon exchange 

c^ ' mechanism tends to cancel the 1/Nc term in the factorized amplitude. 



submitted to Physics Letters B 



1. There has been raising interest during a last few years in testing the factor- 
ization approach to nonleptonic heavy meson decays. The current activity in this 
direction has been triggered by the observation |1|] that the available data for D- 
meson decays seems to support a rule of discarding 1/Nc terms in factorized hadronic 
matrix elements 0, which is in contrast to the standard prescription that keeps such 
terms. This approximate cancellation of 1/Nc terms in the D decays has been ex- 
plicitly checked P| within the QCD sum rule approach. Surprisely, the recent data 
on B-meson decays 0] signals that discarding l/N^. terms is not a universal rule for 
all channels. The l/A'"c-suppression has rather a dynamical character and varies for 
different channels [5-10]. Still, there exist (both perturbative and non-perturbative) 
indications that patterns of deviation from the factorization approximation are alike 
inside each separate class of decays if they are classified according to factorization 
properties of corresponding effective Hamiltonians. In particular, it has been found 
in Ref. [§],|01 that the 1/Nc rule is likely to hold in the Class I decays (see (5) below). 
The effect has a dynamical origin, and is due to non-perturbative gluon effects which 
have been estimated in 0,[^,0 by the QCD sum rule method. The present letter is 
aimed to test an importance of these effects on the so-called color suppressed (Class 
II) decays. As a particular example, we will study the decay B^ -^ D^tt^ following 
the method suggested in Ref. 0. 

At the tree level, weak hadronic Cabibbo favored B-decays correspond to the 
quark transitions b -^ ccs and b -^ cud and are governed by the effective Hamiltonian 

H,„ = ^[Vc,V:,{C,{fi)01 + CMO^) + Vc,V:,{C,{fi)0- + CMO-)] (1) 

where ( T^ = 7^(1 - 75) ) 

or = {cT^b){dT^u) and O^ = {dT^b){cT^u) (2) 

and the operators OJ2 are obtained from 0^2 by the substitution {d,u) -^ (s, c). 
The Wilson coefficients Ci{fi) are due to the renormalization of the bare Hamiltonian 
Hw ~ Oi by hard gluons with virtualities larger than /x^ = O(m^). In the leading- 
log approximation Ci^2 = |(C'+ i C*-) with |ll| 



67± 

with 7_ = -27+ = 2 . At ;u ^ 5 GeV and Ufi = 5, this yields ( Ams ^ 200 MeV ) 

Ci = 1.117 and C2 = -0.266 (4) 

Within the factorization approximation one can distinguish between three classes of 
decays for which the corresponding amplitudes have the following structure 

Aj ~ ai{fi) < Oi >, Ajj ~ 02 (/i) < O2 >, Aijj ~ [ai(/i) + xa2(/i)] < Oi > (5) 



Here < Oi > are the (factorized) hadronic matrix elements of the operators Oi and 
ai{fi) are QCD factors related to the coefficients Cj(yu) : 

ai(/i) = Ciifi) + ^C'2(/i) , a2{fi) = C2{fi) + ^Ciifi) (6) 



An attempt of the global fit of nonleptonic B-decays yields [|12| , |T3 

ai = 1.05 ± 0.10 and a^ = 0.25 ± 0.05 (7) 

It should be mentioned that the very possibility of the global fit is rather questionable 
in view of both experimental and theoretical uncertainties and the expected variation 
between the naive factorization and the 1/Nc rule for different channels [5-10]. Still, 
note that while the value of ai is consistent with suggestive discarding the 1/Nc 
correction in (6), this is not the case for a2. The problem of the (non-)factorization 
thus becomes particularly important for the decays falling into the second class in 
the above scheme. These are called color suppressed decays in view of the fact 
that they are non-zero only due to the QCD-induced operator O2. To illustrate 
the basic difference between the class I and class II amplitudes, consider the decays 
B^ -^ D^Ti^ and B^ -^ D^n^. For the former decay, the operator Oi gives the 
leading contribution a.s Nc ^ 00 : 

(D+7r-|0i|5°) = {n\dT^u\0){D+\cT^b\B) + 0{^) , (8) 

c 

while O2 is factorized only after applying the Fiertz transformation 

O2 = ^Oi + 2{cT,t%){dT^t''u) = ^Oi + O2 (9) 

Note that the matrix element {D^it^\02\B) vanishes in the factorization approxima- 
tion because of the color conservation. An appealing method to estimate this non- 
factorizable contribution has been proposed in 0. There it has been shown that the 
above matrix element can be reduced to a simpler one (p^J^^ — p\^^){D\cV ^gG a^h\B) 
by virtue of the short distance operator product expansion (OPE) , while the latter 
matrix element is fixed by the heavy fiavor symmetry and can be expressed via the 
experimental number M^^ — M^ ~ 0.46 GeV^. Consequently, the non-factorizable 
piece is of the same order as the factorizable (1/A'"c)0i part of O2, and has the 
opposite sign. This observation justifies the (approximate) rule of discarding 1/A^c 
corrections in the class I decays. 

Let us now consider the B^ -^ D^-k^ decay. Then O2 factorizes and 

(DV|ife//|5°) ~ (C2 + ^){D\cV,u\Q) {T:\dV, h\B) 

+Ci{D-K\2{cV^t''u){dV^t''b)\B) (10) 

Our task is to estimate the non-factorizable contribution to this decay which is given 
by the second matrix element in (10). Note that within the approach of Ref. [Q, this 



contribution will be proportional to the matrix element pQ,(7r| JF^Ga^fel-B) (see below) 
, which is not known from general principles and deserves a separate calculation. We 
will present below an estimate for {D^n^\Oi\B^) based on the so-called light cone 
QCD sum rules method ( see |T^ , [|1^ and references therein) which combines the 



traditional technique of the QCD sum rules |T^ with a description of an emitted 
light particle ( 7r° in our case ) in terms of the wave functions of increasing twists. 
We will argue that the suppression of 1/Nc contributions via the soft gluon exchange 
mechanism ^ holds also for 5° — > D^n^ and probably for other Class II decays. 
However, for the color suppressed decays, such cancellation is not what is welcomed 
by the phenomenology if one assumes the validity of the global fit (7). We will come 
back to this point later on. 

2. To implement the two-step strategy suggested in Ref. ||^, we start with the 
correlation function 

A^ = I /rfxe-^-'(7r°|^(x)7^75c(a:),Qi(0)|^°) = ifDP,{D'K\d^\B) , ^ ^ + ... (11) 

J Til j~\ P 

where the ellipses stand for higher states contributions. At large Euclidean p^ the 
correlation function (11) can be calculated in QCD. The leading contribution is due 
to the soft gluon emission from the quark loop. A simple calculation yields the 
following "sum rule" M 



ifj,{DiT\d,\B)—^ = —^p^{rr\dgG,,j,j,b\B) 
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mi) — p^ Att^ 

(12) 
which is to be satisfied in the duality interval 1 GeV^ < —p^ < 4 GeV^. (In 
obtaining (12), we have omitted a contribution proportional to the matrix element 
p^{Tx\dgG ^u'~^yh\B) ~ pQ{T[\dgGQi'~^ih\B) as it is down by the inverse heavy quark 
mass : 7i6 = 0(— ) ). As the ratio of the kinematic factors remains approximately 



•■mf,' 



constant and equal 1 in the duality region, we obtain 

{D\p)T:\q)\d,\B\p + q)) ^ --l^p^{n{q)\dgG,.^.^,b\B{p + q)) (13) 

(Note that in the case at hand the use of the short distance OPE is justified by 
the fact that the final c-quark is heavy and its velocity is small in the rest frame 
of the b-quark |0 ). It is important to point out that the formula (13) is obtained 
at the zero order in a^, and cannot be considered as the renormalization group 
covariant one. The answer (13) refers to the normalization point n = 0{mb), the 
entire /^-dependence being implicit. At the one-loop level the operator Oi mixes with 
O2, while the operator dgG^i^'y^'y^b mixes with the operator m'ld'-y^h and operators 
vanishing on the equations of motion : 



+ -f^\nqi--diWl, + l,V)b (14) 

z 47r Hi o 

where we have set nia = and Cp = ^n'^ ~ ^Z^- Thus, one has to understand 
(13) as the tree level relation. For consistency, in what follows we will omit one-loop 
contributions altogether. 

To find the new matrix element (13), consider another correlation function 

Tc(p,q) = i I dxe"'{7r(q)\d(x)gG,„'f^',sb(x)b(0)i',!,d(0)\0) 

= ^;;43^b'"/.(^V(P" + 2,J/,(p'0l (15) 

where {B\bij5d\0) = m'^fB/mb . At this stage the light cone QCD sum rules method 
suggests a simple and straightforward way of the calculation. At large Euclidean 
{p + qY the leading contribution to (15) is 

/(i ifi h* ~ 

(27r)4(m^ - p^ e-^^"'^^'(vr|J(x)G.^(x)7^75(fc + m,)75C^(0)|0) (16) 

(here k = k^'jfj), that can be further evaluated introducing the pion wave function 
(WF) of twist 3 (p-sn 

— % 

(7r|d(x)5fG'^^(t;x)o-c,/375rf(0)|0) = — y=f37r[qaiqi,guf3 - quQ^ip) - Qf^igt^gua - quQ^ia)] 

X /Z}a,03.(a.)e*^"("^+^"^)(17) 



here Dai = daida2da^5{ai + a2 + Ci^ — l) and f^T, ~ 0.0035 GeV^ for /i^ ^ 1 GeV^ 
13 ) , and the set of WF's of twist 4 : 



{Tx^\d{x)-i^igGap{vx)d{Q)\Qi) = — y= 



qpiQa^^ ) - qa{9|3^, , 

qx qx 



X f Dai4)±{aie'''^^'''+'''''^ 



-Ij=^{q^xp - qpx^) /Z}«,0||(a,)e^''^("^+^"^) (18) 
Y / qx J 

Two more WF's of twist 4 0^ , 0|| are defined analogously to (17) with the substi- 
tution (igGap) -^ {l5gGaf3) )• 



Then, using the identity 

/ du Dai5{u — ai — as)^{ai) = du da3^{ai = u — a^, a2 = u, a^) , (19) 
Jo J Jo Jo 

the answer for the simplest tree diagram can be written as 

^" = ~^ / ~~2 7 \ ^ / C?"3[-2(pg)/37r037r + fnmb{<l)\\ - 20_l)](m " 03,^,03) 

(20) 
A systematic study of higher twist WF's beyond the asymptotic regime has been 
done in Ref . ||19[ . This analysis has been based on the expansion in representations 
of the so-called collinear conformal group S0(2,l), which is a subgroup of the full 
conformal group SO (4,2) acting on the light cone. The asymptotic WF's are defined 
as contributions of operators with the lowest conformal spin and unambigouosly fixed 
by the group structure. Pre-asymptotic corrections correspond to the operators with 
the next-to-leading conformal spin, whose numerical values have been calculated by 
the QCD sum rules method. The result for <^3^ reads [^ : 

037r(aj) = 360a;ia;2a3[l + cJi,o^(7a3 - 3) -|-u;2,o(2 - 4aia2 - 80:3 + Sag) 

+ uji^i{3aia2-2a3 + 3al) + ...] (21) 



where tuifi = —2.88 , 1^72,0 = 10.5 , uji^i = in a low normalization point [0,||T9|. 
For the twist 4 WF's the relevant first order formulas are 

0x(«,) = 30(52(1 -«3)«2[1^2£(l-2«3)] 

4>\\{ai) = -1206^a^a2a3[^ + e{l-3a3)] (22) 

where the parameter 5"^ ^ 0.2 GeV^ ( at yU^ ~ 1 GeV^ ) is defined via 

{0\dgG^u7u75u\'^^) = iS^fnqf, (23) 

and e ~ 0.5 is the weight of the first conformal spin correction. We will need the 
above set of parameters renormalized to a higher normalization point /x^ ~ /i^ = 

771^ — ml ~ (2.4 GeVy which is the characteristic virtuality of the b-quark in 



the B-meson 1T3],[|15|. The corresponding anomalous dimensions can be found in 
ig. The resuhs read ^ f^^ = 0.0026 GeV^, wi,o = -2.18, cui,i = -2.59, uj2,o = 
8.12, 52-0.18, £-0.4. 

To match the answer (20) with the B-meson contribution to the correlation 
function (15), we note that (20) can be re- written as the dispersion integral with 
the expression {ml — up'^)/u being the mass of the intermediate state. The duality 
prescription tells that this invariant mass has to be restricted from above by the 
duality threshold Sq — 35 GeV^ (this value is obtained from corresponding two- 
point sum rules). As it is easy to see, this transforms into the effective cut-off from 



below in the u-integral [^|,|T^. Finally, we make the standard Borel transformation 
suppressing both higher states resonances and higher Fock states in the full pion 
wave function. Under the Borel transformation —{p + q)"^—^ M"^ 



m 



m. 



{p + qy 



exp 



B- 



M2' 



1 



ml — {p + uqY 
Our final sum rule takes the form 

1 



1 2—2 

1 , mf — up , 



(24) 



h{p' 



-f2{p' 



nib 



^ n , I du da'i exp — -— — - — B m 



-JStt 2 '^Svr + In (0|| - 



(tti) 



(25) 



3. Now we turn to numerical estimates. In evaluating (25), we have used the 
following set of parameters : nib = 4.7 GeV, m^ = 5.28 GeV, Sq ~ 35 GeV"^, fs — 
135 MeV ||2^,|]l5l. The Borel mass M^ has been varied in the interval from 8 to 
20 GeV^. We have found that within the variation of M^ in this region, the result 
changes no more that by 10 % and yields for p^ '. 



^l 



fi{ml) ~ -/2(m^) 



1 

V2 



X 0.08 GeV^ 



Then for the matrix element of interest we obtain 
pMdGaf.l,l5b\B) ~ Mlh{ml) + MlUml) 



V2 



X 1.9 GeV^ 



(26) 



(27) 



For the factorizable amplitude due to the operator (l/A^c)C'2 we have 

1 

i 

Nr. 



M 



f 



-^fDP^,{^r{q)\dJ^b\B{p + q)) 
-fDP,[2q,fHp')+p,{fHp') + /."(/))] 



(2^ 



The value of the form factor f^irrvj-,) can be read off the results of Ref. ||20| where 
this quantity has been calculated (for the charged pion) by virtue of the light cone 
QCD sum rules method. The answer is 



f^rnl) 



'V2 



X 0.3 



(29) 



where the number (— 1/-\/2) is due to the different isospin structure in our case. For 
the second form factor we use the model [01 



-^„2\ 



f;ip 



-UiP - -fniP 

niB + ni^ 



(30) 



and therefore neglect the second term in (28). In the nomenclature of Ref. [Q, 
we thus obtain the following estimate for the ratio of the non-factorizable to the 
factorizable 1/iVc amplitudes : 

- Mf - A7^^fl{ml-ml)ft{ml)- "'^ ^' ^ 

where we have used the value //) ~ 170 MeV corresponding to omitted a^-corrections 



in the relevant two-point sum rule [g^; [HI- Thus, our final result (31) suggests that 



the non-factorizable contribution tends to cancel the factorizable 1/Nc amplitude 
due to the operator {1/Nc)02 (see (10)), in agreement with expectations of Ref. 
|^,[^. The sign of the effect can also be compared with the estimate within the 
QCD sum rules method done in Ref. for the weak decay B -^ J/"^ K, which 
also belongs to the class of color suppressed decays. The authors of have found 
that the power corrections due to the gluon condensate partly cancel the 1/Nc fac- 
torizable amplitude, i.e. their sign of r is negative, too. On the other hand, within 
the approach of Ref. 0] the non-factorizable amplitude for this decay can be ap- 
proximately expressed via the matrix element pa{K\sGafi'jfi'y5b\B) , which can be 
extracted from (27) assuming the SU(3) limit. One has to bear in mind, however, 
that the approach of Ref. cannot be literally applied to the decay B —>■ J /'^ K. 
The poor stability of the corresponding sum rule |]^ implies that perturbative cor- 
rections or operators of higher dimensions must be important there. 

4. In this letter we have estimated the soft gluon exchange mechanism contri- 
bution to the non-factorizable amplitude of the B^ — > D^vr" decay and found the 
tendency for the cancellation of 1/iVc, in seeming contradiction with (7). One should 
emphasize, however, that our result does not mean the actual contradiction of the 
theory with the available data [^. One possible sourse of the disagreement may 
be the unjustified use of the global fit leading to the particular numbers (7). Such 
fit implies that the factorization properties are alike in all non-leptonic two-body 
B-decays. The validity of this assumption has to be examined in QCD. Probably 
the more important origin of disagreement with (7) are large perturbative 0{as)- 
corrections which are not taken into account in our calculation, and correspond 
to different contributions to non-factorizable amplitudes. Estimates made in Ref. 
indicate that radiative corrections are important for a complete evaluation of 
non-factorizable amplitudes for the Class II decays. As we have mentioned, in a cal- 
culation including radiative corrections one needs an accurate separation of 0{as) 
corrections related to the matrix element itself and 0{as) terms due to the mixing 
with the two-particle operators, cf. (14). 

In the recent paper [|I^] it has been stressed that non-factorizable amplitudes 
are very important in one more aspect. The point is that the coefficient a2 becomes 
strongly /x- and scheme- dependent beyond the leading-log approximation. Depend- 
ing on the renormalization scheme, 02 can scale crudely from 0.1 to 0.2. This is in 
strong contrast with oi which exhibits very weak /z— and scheme- dependence. The 
importance of higher order QCD corrections for an accurate calculation of 02 can 



be intuitively understood as a consequence of the fact that 02 is the difference of 
the nearly equal numbers, and thus is very sensetive to their precise values. Since 
the factorized amplitudes are /i— and scheme- independent, only the non-factorizable 
contributions can remove this scheme (and /i- ) dependence in the physical ampli- 
tudes. Unfortunately, we have nothing to add in respect to this important problem 
in view of our neglecting perturbative 0{as) corrections. At the same time, the 
observation of Ref. 0] suggests that the perturbative corrections to the matrix 
elements are presumably more important for the Class II decays than for those of 
the Class I. The main question is whether the account for hard gluon loops is able to 
change the sign of the ratio (31). A calculation of the matrix element (D^ir^lOilB^) 
including the radiative corrections can be done either by the methods of Ref. 0, 
0, 0, or directly by the light cone QCD sum rules method |T^. In the latter 



approach, one has to calculate the three-point correlation function of the D- and 
B-meson currents and the effective Hamiltonian between the vacuum and the pion 
states. We have explicitly checked that in this case the leading 0{a^) contribution 
is again given by a combination of the three-particle WF's of twists 3 and 4. How- 
ever, the corresponding estimate based on retaining only these terms manifests a 
poor stability, that indicates the importance of radiative corrections or/and higher 
twist effects in the corresponding sum rule. This instability does not occur in the 
above sum rule, that formally justifies our neglecting the radiative and higher twist 
correction. Undoubtfully, the complete evaluation of perturbative gluon effects on 
hadronic matrix elements of interest is needed before any comparison with the data. 
To our knowledge, this work is currently in progress 0. 

I wish to thank Boris Blok for drawing my attention to this problem and valuable 
discussions. 
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